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This report summarizes our study of Neutral Current (NC)-induced photon production in Mini- 
BooNE, as motivated by the low energy excess in this experiment [A. A. Aquilar-Arevalo et al. 
(MiniBooNE Collaboration), Phys. Rev. Lett. 98, 231801 (2007); 103, 111801 (2009)]. It was 
proposed that NC photon production with two anomalous photon- Z boson-vector meson couplings 
might explain the excess. However, our computed event numbers in both neutrino and antineutrino 
runs are consistent with the previous MiniBooNE estimate that is based on their pion production 
measurement. Various nuclear effects discussed in our previous works, including nucleon Fermi 
motion, Pauli blocking, and the A resonance broadening in the nucleus, are taken into account. 
Uncertainty due to the two anomalous terms and nuclear effects are studied in a conservative way. 

PACS numbers: 25.30.Pt; 12.15.Ji; 14.60.Lm 
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I. INTRODUCTION 



The MiniBooNE neutrino experiment searches for 
Vfj, v e and — > v e oscillations as motivated by 
the anomalous excess of v e observed by LSND in the Sy, 
beam Q. Such oscillations were not clearly observed in 
MiniBooNE's neutrino run, but seem to exist in the an- 
tineutrino run, hinting at possible CP violation. At the 
same time, MiniBooNE reported a low energy excess in 
both runs |3|. It is natural to ask whether the two are 
interrelated. In other words, it is important to under- 
stand this excess in a complete neutrino oscillation anal- 
ysis. For the excess, different explanations have been hy- 
pothesized, including anomalous photon- ,Z boson- vector 
meson couplings Lorentz violation Q, sterile neu- 

trinos and new gauge bosons @. However, a clear 
understanding of the excess has not been achieved so 
far. It was pointed out in Refs. [1] that neutral current 
(NC)-induced photon and 7r° production are potential 
backgrounds in the experiment, because the final photons 
can be misidentified as electrons in the detector. Mea- 
surement of neutrino-induced pion production has been 
carried out in MiniBooNE and used in their background 
analysis but NC photon production has not been 
directly measured. Instead, NC photon production has 
been estimated by assuming a relationship between A 
radiative decay and pion decay in the nucleus (Tlj . On 
the other hand, the authors in Refs. 4-6] proposed that 
anomalous photon- Z-meson couplings might increase NC 
photon production significantly to explain the excess. 



A rigorous calculation of NC photon production, tak- 
ing into account nuclear effects, has not been carried out 
to date but is what we aim to accomplish in this study. 
There arc different channels in photon production, in- 
cluding incoherent and coherent production, but such 
distinction is not obvious in MiniBooNE's background 
analysis. The nucleon Fermi motion, Pauli blocking, and 
broadening of the A inside nucleus [due to the opening of 
other decay channels (e.g. AA" — > NN)] have important 
effects on the production, but unfortunately are missing 
in the previous calculations This report summa- 

rizes our computation of NC photon event numbers in 
MiniBooNE. The calculation is based on a series of works 
dedicated to low energy (neutrino energy below 0.5 GeV) 
neutrino-induced photon and pion production fl2l-[l5] in 
which the neutrino-nucleon interaction kernel as derived 
from a chiral effective field theory (EFT) for nuclei 1 is 
calibrated to available pion data. Also, the nuclear effects 
and the approximation schemes in the study of incoherent 
and coherent production are benchmarked by the corre- 
sponding electron-nucleus scattering and pion photopro- 
duction data. Two contact terms in the EFT, named as 
Ci and ei terms here, are the low energy manifestation 
of the prop osed anomalous photon- Z-meson interactions 
[1 03, 03]- Furthermore, it was realized Q3, 03] that 
other sources also lead to these low energy interactions. 
In Ref. [HI, we have calibrated the ci and e\ values to 
the coherent pion photoproduction data, and found C\ 
is correlated with the other two parameters (r s , r v ) that 
control the A binding potential and spin-orbit coupling in 
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The EFT is knowns as quantum hadrodynamics EFT (QHD 
EFT). The motivation for this EFT and some calculated results 
are discussed in Refs. |16H23H . 



the nucleus. The uncertainty of these parameters, which 
eventually leads to uncertainty in the photon event cal- 
culation, will be discussed in detail later. 

As emphasized in the previous works [l2Hl"5| , the EFT 
approach is only valid for low energy neutrinos. It turns 
out that the so-called low energy excess, if due to NC pho- 
ton production, gets important contribution from higher 
energy neutrinos because the reconstructed neutrino en- 
ergy Eqe in MiniBooNE, based on quasielastic scattering 
kinematics, can underestimate the true energy substan- 
tially in this channel. This study essentially amounts to 
the high energy extrapolation of previous works [1 2Hl 51 ) . 
To benchmark the extrapolation of the neutrino-nucleon 
interaction kernel, the neutrino-induced pion production 
is calculated and compared to higher energy neutrino 
data (~ 1 GeV). A reasonable form factor is associated 
with ci and e\ couplings to regularize their high energy 
behavior as inspired by the anomalous photon- Z-meson 
interactions [||. The nuclear effects are included as be- 
fore. The final results, based on the detector information 
and updated detection efficiencies from Refs. 0, |2fj, in- 
dicate that in both neutrino and antineutrino runs, for 
Eqe below 0.475 GeV, the experimental estimate is close 
to our computation; at higher Eqe, it is close to our lower 
bound but significantly smaller than our upper bound, 
although still compatible within the experimental un- 
certainty. We also compare our event numbers vs. E y 
with the experimental estimate: our results agree with 
the experimental estimates, except in the lowest E n bin 
([0.1 , 0.2] GeV) in both runs where ours are larger than 
their estimates yet the differences are not large enough to 
explain all the excess. In other words, the excess needs 
to be resolved from other perspectives. Furthermore, the 
increase in the number of events due to c\ and e\ at 
higher Eqe is more substantial in antineutrino than in 
neutrino data, which raises another interesting point: the 
neutrino and antineutrino respond differently to the two 
couplings. It is possible that, similar to photon produc- 
tion, the mechanism responsible for the excess can affect 
the event count differently in the two runs. Hence the 
statement [l| should be taken cautiously that the data in 
the antineutrino run above Eqe = 0.475 GeV is with- 
out any unknown background solely because the neutrino 
run indicates so. 



The rest of this report is organized as follows. In 
Sec. HH we show our high energy extrapolation of 
neutrino-induced pion production and compare it with 
available data. After that, Sec. IIIII summarizes our com- 
puted photon event number distributions against differ- 
ent kinematic variables, and comparisons are made be- 
tween them and the MiniBooNE estimates. The uncer- 
tainty of our computation is discussed in detail there. 
Finally in Sec. IIV1 the report ends with a short summary 
and discussion. 




(a) (b) (c) (d) (e) (f) 



FIG. 1: The Feynman diagrams for photon production in- 
duced by various currents shown as "C" in the plot including 
vector, axial, and baryond currents. The diagrams for pion 
production can be viewed as these diagrams with final pho- 
ton lines changed to pion lines. Diagrams (a) and (b) have 
the A as intermediate states; (c) and (d) have the nucleon as 
intermediate states; (e) has a pion pole in the t channel; (f) 
includes the contact terms. 



II. HIGH ENERGY EXTRAPOLATIONS AND 
BENCHMARKS 

First, to benchmark our formalism, we apply it to the 
study of neutrino-induced pion production from free nu- 
cleons with neutrino energy up to ~ 2 GeV, and check 
it against the existing data from ANL and BNL exper- 
iments [27], The contributing Feynman diagrams 
are similar to those in Fig. [T] with the final photon lines 
changed to pion lines. See Ref. [HI, [l3| for a detailed dis- 
cussion of these Feynman diagrams included in our previ- 
ous low energy pion production study. However, only the 
contributions from the A resonance and nucleon interme- 
diate states [diagrams (a)(b)(c) and (d)] are considered 
here. Several contact terms (besides c\ and ei) exist, 
whose high energy behavior is unknown, and hence are 
not discussed here. To extrapolate our previous calcula- 
tions [13I - H5I ] to higher energy, we make use of the fitted 
form factors of the axial transition currents (N O A) 
in Ref. [3l| . The form factors of the vector transition 
currents are constrained by the corresponding Electro- 
magnetic interactions 0, [3(| • The form factors of the 
nucleon currents [as used in diagrams (c) and (d)] are 
from Ref. [32| (for the vector current) and Ref. (33| (for 
the axial current). Fig.[2]plots charged current (CC) pion 
production data from the ANL [27j and BNL [28| exper- 
iments. The former requires the final pion and nucleon 
center of mass energy M nn to be less than 1.4 GeV, but 
no such cut is applied in the BNL data. Two different 
calculations, with and without such a cut, are shown. 
We see that the "no cut" computation agrees reasonably 
well with the BNL data, while the "with cut" slightly 
overestimates the cross sections as indicated by the ANL 
data, which strongly indicates that the contact terms if 
regularized cannot be a dominant contribution in the ~ 1 
GeV region. 

Second, let's proceed to NC photon production from 
nucleons by using the mentioned form factors and the 
procedure handling them [l3|. Fig. Q] shows a list of 
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FIG. 2: Total cross section vs. neutrino energy in CC pion 
production from free nucleons. In the ANL data [27j | . the final 
pion and nucleon center of mass energy M nn is required to be 
less than 1.4 GeV, while no such cut is applied in the BNL 
data [2^1 . Two different calculations are shown corresponding 
to with and without the cut. 



the diagrams. However, we can neglect the contributions 
from diagram (e) because of the 1 — 4 sin 2 6 W suppression 
(it contributes zero in coherent production). In diagram 
(f), two couplings, the ci and e\ terms, contribute here. 
It was pointed out in Ref. [5| that the two terms can be 
induced by the w and p anomalous couplings with the 
photon and Z boson (ci ~ 1.5 and e\ ~ 0.8), but also 
emphasized in Ref. [EEl that other sources exist, includ- 
ing A off-shell interactions. At low energy, we calibrated 
ci to be around 1.5 ~ 3 by using coherent pion photopro- 



duction data [T5|. The two terms contribute little in the 
production from nucleons and the incoherent production 
with neutrino energy up to 0.5 GeV [El [T3I fl3 . but can 
be relevant in the coherent production [5j, Il5| . To reg- 
ularize their high energy behavior, we use the following 
form factors j5|: 



(Pf-Pi) 2 



(1.0 GeV) 2 

where pf — pi is the momentum difference of the scattered 
nucleon in the interaction kernel. The total cross sec- 
tions for the production from free protons and neutrons 
are shown in Fig. |3l For example, "p(f)" corresponds to 
the production from free protons, while "p(b)" is for the 
incoherent production from bound protons in 12 C, which 
will be addressed later. In each case, three calculations 
including different diagrams are shown: "only A" has di- 
agram (a) in Fig. [1} "A + N" includes (a) (b) (c) and 
(d); "full" includes all the diagrams [(e) is neglected as 
explained above]. From now on, we focus on the ~ 1 
GeV region because (anti)neutrinos with higher energy 
are suppressed in the MiniBooNE spectrum. Clearly, the 
A dominates in both neutrino and antineutrino scatter- 
ings. For the proton, adding N contributions increases 
the cross section by ~ 10%, while for the neutron its 
contribution is less. In the "full" calculation, adding the 
contact terms, c\ — 3.0 and e\ = 0.8, does not change 
the results at low energy ~ 0.5 GeV, but increases the 
cross section well above ~ 1 GeV (they are more signfi- 
cant in the antineutrino-induced than neutrino-induced 
productions). Since the contact terms and the associated 
form factors are not well understood so far, the increases 
well beyond ~ 1 GeV should be taken cautiously. 

Third, to calculate the productions in neutrino-nucleus 
scattering, the medium-modification of the interaction 
kernel needs to be included. The ground state of the 
nucleus (e.g. 12 C) is computed in the mean-field ap- 
proximation of QHD EFT [II HI]. We apply the Local 
Fermi Gas approximation to calculate incoherent produc- 
tion and the "optimal approximation" for coherent pro- 
duction [l4l [THj ] . The medium modifications have been 
checked against the inclusive electron- 12 C scattering and 
coherent photoproduction data [II],[l5j]. In principle, the 
available data [25[ for (in) coherent neutrino-induced pion 
production from nuclei can be used to constrain theoret- 
ical models. However, before doing so, final state inter- 
action of pions and nucleons have to be studied carefully, 
which is not the focus here. In addition, for coherent 
pion production, several contact terms besides c\ [(f) in 
Fig. [1] with final photon line changed to pion line] can be 
important Il5ll. but their high energy behavior are not 
clear so far |l3tl34|. Fortunately, for photon production, 
the final state interaction can be ignored, and only two 
contact terms exists (as motivated by the anomalous cou- 
plings) , which makes the high energy extrapolation more 
tractable for neutrino energy up to 1 GeV. The total 
cross sections for incoherent NC photon production are 
shown in Fig. [3] in comparison with those from free nu- 
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cleons. The labeling of the different curves have already 
been mentioned above. Two parameters, (r s , r v ), which 
are the meson-A couplings and related to the A bind- 
ing and spin-orbit coupling, are set to be (1, 1) in these 
calculations. The incoherent production is not sensitive 
to different r s and r v values around (f , f) (see Ref. [lH 
for the details). Clearly, the A contribution is reduced 
by 50% compared to that from free nucleons in all the 
different channels. Because of the opening of other decay 
channels (e.g. AN -s- NN), the A width increases in the 
nucleus. As a result, the flux in the radiate decay channel 
is decreased. Meanwhile, the nonresonant contributions 
from both N and contact terms are reduced less, so the 
total cross sections for the full calculations are reduced 
less (~ 30%). 

In Fig.EJ we show the calculated total cross sections for 
coherent photon production in (anti)neutrino- 12 C scat- 
tering. "(1, I)" for example indicates (r s , r v ) = (1, 1) 
as mentioned before. Results of using two different sets 
of diagrams are compared: A + N includes contributions 
from both A and N, while the full calculation has all 
the diagrams. We can see the coherent production is in- 
deed sensitive to (r s , r v ) and c\. The choices of them 
in the two "full" calculations are due to the calibration 
in Ref. based on coherent pion photoproduction (ei 
does not contribute in the isoscalar nucleus). Interest- 
ingly, for each calculation the difference between the neu- 
trino and antineutrino-induced production cross section 
is reduced with increasing (anti)neutrino energy, which 
is consistent with previous discussion [15| . 



III. EVENT NUMBERS 

In this section, we summarize our computation of the 
NC photon event numbers in MiniBooNE. The spectrum 
of Vfj, and can be found in Refs. [24| : the median energy 
is around 0.5 ~ 1.0 GeV, and the maximum energy is up 
to ~ 2.0 GeV. The proton on target (POT) numbers for 
the neutrino (antineutrino) run is 6.46 (11.27) xlO 20 . The 
total effective mass of the mineral oil (CH^) detector is 
8.12 x 10 8 g. The photon detection efficiency vs. photon 
energy is listed in Tab. |H For photon energy higher than 
1.0 GeV, we assume the efficiency to be 0.102, which 
should not cause mistakes because our calculations show 
that most photon events have energy lower than 0.8 GeV. 

Tabs, im and IIIII list the photon event numbers against 
Eqe (in three bins) for neutrino and antineutrino runs 
(the formula for reconstructing Eqe can be found in 
Ref. 24]). The events are broken into different chan- 
nels including coherent ("coh") and incoherent ("inc") 
production, and the production from unbound protons 
("H"). "Total" corresponds to the total number, to which 
MiniBooNE's estimate and excess should be compared. 
In each entry, we list two numbers showing the lower and 
upper bound of our computation. The lower bound in all 
three channels is based on the "A + N" calculation with 
(r s , r v ) = (1, 1), while the upper bound corresponds to 
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FIG. 3: Total cross section (per nucleon) for the photon pro- 
duction from free nucleons and the incoherent production 
from 12 C. The first (last) two are for neutrino (antineutrino)- 
induced production. See the text for detailed discussion of the 
different curves. 
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FIG. 4: Total cross section (per nuclei) for the coherent pho- 
ton production in (anti)neutrino- 12 C scattering. The first 
(last) is for neutrino (antineutrino)-induced production. See 
the text for detailed discussion of the different curves. 



the "full calculation" with c% — 3.0, e\ — 0.8, (r s , r v ) = 
(1, 1) [for the production from free protons, (r s , r v ) are 
not relevant]. See Figs. [3] and S] and the text around 
them for a detailed discussion on the different calcula- 
tions. Comparing different calculations in these plots, we 
can see that the choice of the two bounds are reasonable. 
In both runs for the first two Eqe bins, the difference be- 
tween our upper bound and experimental estimate is not 
significant enough to explain the excess. In the final Eqe 
bin, our upper bound is significantly larger than the esti- 
mate; however, this difference should be taken cautiously 
because our simple regularization of contact terms that 
are mostly responsible for the upper-lower difference have 
not been well benchmarked. Still, considering the large 
uncertainty of the excess in the final bin, our results are 
consistent with the experimental estimate indicating no 
excess in this bin. Another observation, as mentioned in 
the end of Sec. HI is that the contact-term induced event 
increase (in percentage) in the final bin is less in the neu- 
trino run than in the antineutrino run, which means the 
two runs respond to the contact terms in different ways. 

Tabs. IIVI and [V] show the same event numbers as in 
Tabs. HI1 and IIIII but against a different variable, photon 
energy E 7 , in five different bins. The labeling has been 
mentioned above. We can see that the experimental esti- 



£ 7 (GeV) 


0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.9 


Efficiency 


0.089 


0.135 


0.139 


0.131 


0.123 


0.116 


0.106 


0.102 



TABLE I: MiniBooNE's photon detection efficiency vs. pho- 
ton energy [l[ . The listed E 1 value is the center of each energy 
bin. The expermental systematic uncertainty is 15%. 



E QE (GeV) 


[0.2 , 0.3] 


[0.3 , 0.475] 


[0.475 , 1.25] 


coh 


1.5 (2.9) 


6.0 (9.2) 


2.1 (8.0) 


inc 


12.0 (14.1) 


25.5 (31.1) 


12.6 (23.2) 


H 


4.1 (4.4) 


10.6 (11.6) 


4.6 (6.3) 


Total 


17.6 (21.4) 


42.1 (51.9) 


19.3 (37.5) 


MfniBN 


19.5 


47.3 


19.4 


Excess 


42.6 ± 25.3 


82.2 ± 23.3 


21.5 ±34.9 



TABLE II: Eqe distribution of the NC photon events in 
the MiniBooNE neutrino run, comparing our estimate to the 
MiniBooNE estimate [H . 



mates agree reasonably well with our results, except for 
the lowest E 1 bin ([0.1 , 0.2] GeV) in which their esti- 
mates fall below our lower bounds in both neutrino and 
antineutrino runs. This is perhaps due to our inclusion 
of nonresonant diagram contributions, which is not man- 
ifestly included in the experimental estimates. 

Comparing our results with those in Ref. [6|, we can 
see that with the updated photon detection efficiencies 
used here, our results are generally smaller than those in 
Ref. @. In Sec. HI we have shown that nuclear effects 
tend to reduce the total cross section substantially, which 
is also reflected in the event number computation. 



IV. SUMMARY 

In summary, we have calculated the total cross sec- 
tions for NC photon production and computed the pho- 
ton production event numbers vs. reconstructed neu- 
trino energy Eqe and photon energy E 1 for MiniBooNE. 
Two contact terms c\ and e\ (with reasonable form fac- 



E QE (GeV) 


[0.2 , 0.3] 


[0.3 , 0.475] 


[0.475 , 1.25] 


coh 


1.0 (2.2) 


3.1 (5.5) 


0.87 (5.4) 


inc 


4.5 (5.3) 


10.0 (12.2) 


4.0 (10.2) 


H 


1.3 (1.6) 


3.6 (4.3) 


1.1 (2.4) 


Total 


6.8 (9.1) 


16.7 (22.0) 


6.0 (18.0) 


MiniBN 


8.8 


16.9 


6.8 


Excess 


34.6 ± 13.6 


23.5 ± 13.4 


20.2 ±22.8 



TABLE III: Eqe distribution of the NC photon events in the 
MiniBooNE antineutrino run, comparing our estimate to the 
MiniBooNE estimate 
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E-,(GeV) 


coli 


inc 


H 


Total 


MiniBN 


Excess 


[0.1 , 0.21 


0.72 (1.5) 


14.0 (15.0) 


4.4 (4.6) 


19.1 (21.1) 


10.6 


52.5 


[0.2 , 0.3] 


3.2 (5.5) 


22.7 (25.2) 


7.8 (8.5) 


33.7 (39.2) 


32.5 


61.2 


[0.3 , 0.4] 


3.7 (5.4) 


12.7 (15.0) 


5.0 (5.6) 


21.4 (26.0) 


24.7 


58.4 


[0.4, 0.5] 


1.0 (1.7) 


5.4 (7.3) 


2.1 (2.4) 


8.5 (11.4) 


12.7 


-9.7 


[0.5 , 0.6] 


0.32 (1.0) 


2.3 (3.9) 


0.75 (1.0) 


3.4 (5.9) 


4.4 


10.5 



TABLE IV: E y distribution of the NC photon events in the MiniBooNE neutrino run, comparing our estimate to the MiniBooNE 
estimate [lj. 



E 1 (GeV) 


coh 


inc 


H 


Total 


MiniBN 


Excess 


[0.1, 0.2] 


0.55 (1.2) 


4.9 (5.5) 


1.4 (1.6) 


6.9 (8.3) 


4.3 


18.8 


[0.2, 0.3] 


2.0 (3.8) 


8.7 (10.3) 


2.9 (3.3) 


13.6 (17.4) 


14.3 


22.6 


[0.3 , 0.4] 


1.8 (3.0) 


4.0 (5.4) 


1.5 (1.8) 


7.3 (10.2) 


9.1 


11.5 


[0.4, 0.5] 


0.36 (1.0) 


1.3 (2.6) 


0.43 (0.66) 


2.1 (4.3) 


3.6 


18.7 


[0.5, 0.6] 


0.10 (0.72) 


0.51 (1.7) 


0.14 (0.36) 


0.75 (2.8) 


1.1 


8.4 



TABLE V: Ey distribution of the NC photon events in the MiniBooNE antineutrino run, comparing our estimate to the 
MiniBooNE estimate [J]. 



tors) that are partially related to the proposed anoma- 
lous photon- Z-meson couplings are also discussed here. 
The other two parameters (r s , r v ) that control the A 
binding and spin-orbit coupling in the nucleus are an- 
other source of uncertainty in our calculations. Previ- 
ously based on inclusive electron-nucleus scattering data 
and phenomenological fits of the A spin-orbit coupling 
in the study of pion-nucleus scattering, we set (r s , r v ) 
to be around (1, 1). The further study of coherent 
pion photoproduction shows fixing (r s , r v ) is entangled 
with c\\ Two sets, (r s , r v ) = (1, 1), c\ = 3.0 and 
( r s, i°v) — (1; 0.9), cj = 1.5, are allowed. 

To test the neutrino-nucleon interaction kernel with 
neutrino energy ~ 1 GeV, the corresponding pion pro- 
duction with A and N contributions included has been 
shown to be in good agreement with available data. 
Based on the transition and other form factors as bench- 
marked in pion production, we calculate the total cross 
section for photon production, including contributions 
from free nucleons, incoherent production, and coher- 
ent production in the high neutrino energy region. As 
we can see, the regularized c\ and e\ contributions are 
relatively small compared to A + N in all the photon 
production channels until neutrino energy goes beyond 
~ 1 GeV, which indicates our full calculation with the 
two contact terms provides a conservative upper bound 
in the event count. We treat the A + N contribution 
as the lower bound. Based on the two bounds, we con- 
clude that the difference (our count — MiniBooNE's es- 
timate) is not significant enough to explain the excess 
with Eqe below 0.475 GeV in both neutrino and an- 
tineutrino runs. For Eqe above 0.475 GeV, although 
our result (upper bound) is significantly larger than the 
estimate, it is still compatible with no excess within the 



given uncertainty range. For the E 7 distribution, again 
the experimental estimates agree with our results, except 
that in the lowest bins in both runs, ours are larger than 
their estimates but the differences could not explain all 
the excess in those bins. Hence, other perspectives need 
to be explored to understand the low energy excess in 
MiniBooNE. In addition, the two contact terms provide 
an exemplary mechanism to which the neutrino and an- 
tineutrino runs in MiniBooNE respond differently. As 
a result, the nonexistence of low energy excess beyond 
0.475 GeV in the former does not immediately indicate 
it is also true in the latter and vice versa. 

To eliminate the uncertainty in the neutrino-nucleon 
interaction kernel, a model independent approach needs 
to be implemented, for example the dispersion analysis, 
which is currently begin pursued. Moreover, this is also 
relevant to the known 7 Z box diagram in the study of 
parity violation of electron-nucleon (nucleus) scattering. 
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